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We present a comparative seismic study of conditions around and beneath isolated sunspots. Using the European Grid 
of Solar Observations' Solar Feature Catalogue of sunspots derived from SOHO/MDI continuum and magnetogram data, 
1996-2005, we identify a set of isolated sunspots by checking that within a Carrington Rotation there were no other spots 
detected in the vicinity. We then use level-2 tracked MDI Dopplergrams available from SOHO website to investigate 
wave-speed perturbations of such sunspots using time-distance helioseismology. 
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1 Introduction 

Local helioseismology provides a powerful diagnostic of 
the solar interior, allowing us to infer from observations of 
acoustic waves the structures and flows beneath the Sun's 
surface. In this work we focus on investigation of wave- 
speed anomalies underneath several a and j3 sunspots se- 
lected using the European Grid of Solar Observations Fea- 
ture Catalogues (Zharkova et al 2005) . 

Previous estimates of sunspot properties from travel- 
time measurements indicate a negative perturbation to the 
sound speed directly underneath sunspots up to the depths 
2-5 Mm, (Kosovichev et al 2000) . Beneath this depth the 
sound speed perturbation is positive. 

2 Data and method 

We have used MDI level-2 tracked Dopplergramdata. Using 
the EGSO Sunspot Solar Feature catalogue we have identi- 
fied a number of sunspots that have no other sunspots within 
a radius of 10 heliographic degrees. Locating the available 
MDI level-2 data we have selected 6 sunspots comprising 
Active Regions NOAA 7973, NOAA 8038, NOAA 9056, 
NOAA 9138, NOAA 9142, NOAA 9779. From MDI data 
we have extracted 5 12-minute long Dopplergram datacubes 
that were processed using time-distance helioseismology 
methods. In addition, to investigate the solar activity in 
these active regions we use EGSO SEC sunspot data in con- 
junction with level- 1.8 MDI Continuum and Magnetogram 
daily synoptic images and NOAA data available via Solar 
Monitor (http://www.solarmonitor.org). Some basic statis- 
tics available for the regions are presented in Table 1 , while 
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information on evolution of each region can be found in Ta- 
ble 2 (see Section 3 for more details). 



2.1 Time-distance method 

The raw data are measured in the form of Dopplergrams 
representing line-of-sight velocity in the lower solar atmo- 
sphere. The pre-processing stage consists of applying high- 
pass and f-mode filters (Gizon and Birch 2005) to re- 
move supergranulation and f-mode data. Next in order to se- 
lect specific acoustic wave-packets and obtain a clear cross- 
coiTelation signal for p-modes phase-speed filtering is ap- 
plied coiTesponding to a number of skip-distances. Cross- 
coiTelation of solar Dopplergram data at two points over 
time (Duvall ef flZ 1993) produces wave-packet like struc- 
tures, coiTesponding to the propagation between the two 
points of packets formed by acoustic waves inside the Sun. 

To avoid stochastic noise, averaging is necessary for ob- 
taining good cross-correlations. In this work, we proceed 
by computing the cross-correlation of a central point with 
a surrounding annulus. The annulus thickness is taken to be 
around 4.5 Mm. We average individual point-to-point cross- 
correlations C(ri , r2 , t) over an annulus centered on ri and 
with a radius A = |ri — r2|. Thus we estimate a point-to- 
annulus cross-coiTelation function C(r, A, t) (Couvidat et 
al 2005) . 

2.1.1 Travel Times 

From point-to-annulus cross-correlation we derive the one- 
way travel times rj/o(r, A) of wavepackets propagating in- 
ward, i, from the annulus to its center and outward, o, by 
fitting separately the negative and positive time-lag portions 
of the cross-coiTelation with Gabor wavelets (Kosovichev 
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NOAA# 


AR Type 


Date & Time of Observation 


Longitude 


Latitude 


Area 


7973 


a 


1996-06-24T13:05:35.559 


43.526 


9.196 


2.37 


8038 


a 


1997-05-1 1T07:08:35.354 


140.854 


20.633 


3.01 


9056 


Q 


2000-06-26T06:23:32.626 


177.125 


-13.113 


7.140 


9138 


Q 


2000-08-24T05:59:31.318 


117.515 


-30.731 


1.44 


9142 


Q 


2000-08-27T17:35:31.359 


46.371 


15.975 


1.41 


9779 




2002-01-13T06:23:30.128 


269.551 


28.526 


4.31 



Table 1 Sunspots used in this study Usted by NOAA numbers. Date, Carrington longitude, latitude and area are taken 
from the EGSO SFC; times are selected to be closest to the available Dopplergram data. Area is measured in heliographic 
degrees squared. Carrington longitude and latitude are given in degrees. 



and Duvall 1997) . We define the mean travel-times r(r, A) 
as r(r, A) = (r,(r, A) + To(r, A))/2. 

As a first approximation, perturbations 5t{y, A) in these 
mean travel-times are linearly related to the sound-speed 
perturbations 6c in the wave propagation region, 

(5T(r,A)= j J dr' j\zK{r~r',z;A)^{r',z),il) 

where S is the area of the region, d is its depth and ST{r, A) 
is defined as the difference between the measured travel- 
time at a given location r on the solar surface and the av- 
erage of the travel-times in the quiet sun. The sensitivity 
kernel for the relative squared sound speed perturbation is 
given by K. 

In this work we use wave-speed sensitivity kernels esti- 
mated using the Rytov approximation (Jensen and Pijpers 
2003) . The Rytov and ray approximations give very simi- 
lar results, with former being more reliable inverting deeper 
structures (Couvidat et al 2004) . 

2.1.2 Inversion 

The code for travel time inversion is based on a multi- 
channel deconvolution algorithm (MCD) (Jensen et al 
1998) enhanced by the addition of horizontal regulariza- 
tion. Since the sensitivity kernels used here are translation- 
ally invariant, the inversion can be considered as a convo- 
lution product. Thus, to speed up the inversion procedure, 
computations are performed in the Fourier domain. We dis- 
cretise equation O following (Jensen 2001) : 

d, = 6T{k,A,), (2) 

=X(fc,z,;A,), (3) 

nij = Ss{k, Zj). (4) 

For each horizontal wavevector k we find the vector m that 
solves 

mm{\\{d-Gm)\\l + e{kr\\Lm\\l}, (5) 

where L is a regularisation operator In this work we take 
L = diag(S(^]^, . . . , s^^J, where so,i the reciprocal of the 
sound-speed in the i-th layer of the reference model. We 
take e(fc)2 = eg(l + \k\^)P withp = 70 and = 10*. We 
invert for 14 layers in depth located at (.36, 1.17, 2. 11, 3.28, 
4.74, 6.41, 8.6 11.22, 14.28, 17.78, 21.79, 26.31, 31.41, 
36.95) Mm. 



3 Discussion 

The results of sound-speed inversions are presented in Fig- 
ures 1-5. We define the point of maximum absolute flux 
in the MDI magnetogram as the center of the sunspot for 
each of the features considered here. The top two panels in 
Figure 1 show the remapped SOHO/MDI magnetogram and 
continuum images for NOAA 8038 with marked lines pass- 
ing through the center of the sunspot representing longitudi- 
nal and latitudinal sound-speed perturbation depth cuts. The 
lower two panels present the actual depth cuts through these 
lines. Figure 2 contains selected sound-speed perturbation 
maps for some of the other Active Regions considered here. 
Figures 3 and 4 show sound-speed variation at two selected 
depths for each of the regions, with contours of the mag- 
netic lines at 30 and 70 percent of the maximum absolute 
magnetic flux overplotted. Figure 5 shows plots of sound- 
speed variation with depth at the center of each sunspot as 
defined above. 

Table 2 contains extracted results of the inversion and 
feature-tracking for each of the regions under consideration. 
The first column contains NOAA number and Carrington 
Rotation number, the second column presents the history of 
the region, e.g. NOAA 9056 was observed as NOAA 9017 
one rotation earlier, while nothing was observed at the loca- 
tion of NOAA 7973 previously. The last column in the table 
contains additional information regarding further develop- 
ment of the Active Region. For example, NOAA 8038 soon 
after the time of observation develops into a bipolar (3\(3 
Active region, while NOAA 9056 was observed as NOAA 
9096 one rotation later. We note from Tables 1 and 2 that 
the regions we describe as new are observed during the min- 
imum and decaying phases of the solar cycle. 

The results presented in the figures are in line with the 
earlier works by Kosovichev et al (2000), Couvidat et al 
(2005) and others. That is, we see a region of reduced sound 
speed directly underneath the sunspot, with a region of 
positive sound-speed perturbation situated at greater depth. 
From Figures 3 and 4 we can see that these regions appear to 
be situated directly underneath each sunspot for the selected 
depths of 0.36 and 6.41 Mm. 

Analysis of the data in Table 2 suggests that sunspots 
associated with the new Active Regions, i.e. the ones that 
appeared within one Carrington Rotation, show larger am- 
plitude of sound-speed perturbation, Sc/c, than the Active 
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Fig. 1 Latitudinal and longitudinal cuts displayed on MDI magnetogram {top left) and continuum {top right) data for 
NOAA 8038 
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Fig. 2 Sound speed perturbation cut through the center of the sunspot {as lower plots in Fig. 1) for {top, left to right) 
NOAA 7973, NOAA 9056, {bottom, left to right) NOAA 9138 and NOAA 9779. 



Regions that were at least as old as one Carrington rotation 
period. The difference is more pronounced for the regions 
of reduced sound-speed profile. Also, the extent of the re- 
gion of reduced sound speed situated directly underneath 
the sunspot at the point of maximum magnetic flux tends 
to be larger for new sunspots. The validity of such sound 
speed anisotropy below the solar surface may have to be 
further verified by applying an umbra mask test (also known 
as cookie cutter test) which excludes observations taken in 
strongly magnetized regions (Zhao and Kosovichev 2003, 



Hughes et al 2005 and Korzennik 2006). Also, we note that 
the forward problem is currently modeUed based on the ba- 
sic assumption that magnetic field only modifies the effec- 
tive sound speed, while the true picture is likely to be more 
complicated and would include showerglass effect, suppres- 
sion and scattering of p-modes and other effects. 
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NOAA# 


CR# 


AR History 


0-crossing 


min(5c/c) 


max(5c/c) 


Comment 


7973 


1910 


new 


2.47 


-0.282 


0.106 




8038 


1922 


new 


2.11 


-0.219 


0.091 


appearing 


9056 


1964 


NO A A 9017 CR1963 


2.12 


-0.176 


0.078 


#9096CR1965 


9138 


1966 


NOAA9100 CR1965 


2.06 


-0.177 


0.073 




9142 


1966 


NOAA 9105 CR1965 


2.03 


-0.141 


0.077 


# 9201CR1968 


9779 


1985 


new 


2.35 


-0.252 


0.094 





Table 2 Properties of the Active Regions. The zero-crossing depth, i.e. the depth where 6c/ c changes from negative to 
positive, is given in megameters. 



4 Conclusions 

In the present work, using the methods of time-distance 
helioseismology, we have looked at six different isolated 
sunspots at different stages of solar cycle evolution. Using 
European Grid of Solar Observations' Solar Feature Cata- 
logues and NOAA data we have extracted a number of pa- 
rameters and tracked the history of each of the Active Re- 
gions. We have extracted the sound speed perturbation at 
the time of observation as a function of depth directly un- 
derneath each region using Rytov-approximation sensitivity 
kernels. We have used the inversion based on a modified 
Multi-Channel Deconvolution method with regularization 
chosen to be the inverse of the mean sound speed at each 
depth layer. 

Our preliminary results indicate systematic differences 
in subsurface structure of the sunspot depending on the 
stage of its evolution or possibly the phase of the solar cy- 
cle. Further investigations have to be carried out in order to 
verify these results. These include making use of sensitivity 
kernels computed to take into account the specific details of 
our travel-time computations and filtering procedure. Inver- 
sion of travel-time data is an ill-posed problem and further 
improvements to the procedure can be made by including 
a noise covariance matrix and experimenting with different 
choices of regularization operator. 

In addition, further work is to be carried out to verify 
these results by investigating a statistically significant num- 
ber of sunspots. 
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Fig. 5 Sound-speed perturbation, as a function of depth, at the center of the sunspot defined as the point of maximum 
absolute magnetic flux. 
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